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ly-^ , Abstract. Numerical computations have suggested that in causal dynamical 

triangulation models of quantum gravity (CDT) the effective dimension of spacetime in 
the ultra-violet is lower than in the infra-red. In this paper we develop a simple model 
based on previous work on random combs, which share some of the properties of CDT, 
O i' in which this effect can be shown to occur analytically. We construct a definition for 

short and long distance spectral dimensions and show that the random comb models 
exhibit scale dependent spectral dimension defined in this way. We also observe that 
a hierarchy of apparent spectral dimensions may be obtained in the cross-over region 
between UV and IR regimes for suitable choices of the continuum variables. Our main 
result is valid for a wide class of tooth length distributions thereby extending previous 
■ work on random combs by Durhuus et al. 
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1. Introduction 

This work was motivated by observations made in some formulations of quantum gravity 
which we will explain shortly. However, the question as to whether it is possible to define 
consistently a spectral dimension which depends on the scale at which a system is probed 
by a random walk is of more general interest. In this paper we will construct a definition 
for such a spectral dimension and show that there are models which do indeed exhibit 
scale dependent spectral dimension defined in this way. 

The demonstration by t' Hooft and Veltman that quantised General Relativity 
is perturbatively non-renormalisable in four dimensions [T] led to the search for non- 
perturbative formulations of quantum gravity and there are now several approaches 
to this problem. It was first advocated by Weinberg |2J that there might be a non- 
trivial ultraviolet fixed-point and this has been pursued in continuum calculations by a 
number of authors [TO]. [TTj . There is now quite a lot of evidence for such a fixed point 
although it is not conclusive; precisely because GR is perturbatively non-renormalisable 
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it is necessary to make some assumptions or ultimately uncontrolled approximations 
in these calculations. An alternative approach within the fixed-point philosophy is to 
discretise space-time and to look for a critical point or line where a continuum limit 
may be taken to recover continuum gravity. 

Early attempts based on the so-called Euclidean quantum gravity model (see for 
instance [6]) did not lead to a continuum limit in four dimensions but the situation 
improved with the introduction of the Causal Dynamical Triangulation model (CDT) by 
Ambj0rn and Loll in 1998 [4J. This defines the gravitational path integral as a sum over 
discretised geometries (see [5] for a recent review) . In contrast with the earlier Euclidean 
quantum gravity model, the CDT approach takes account of the Lorentzian nature of 
the path integral by building in a well defined temporal structure from the start. As an 
example of the success of the CDT approach numerical simulations have shown that in 
four- dimensional CDT large scale structure in the form of a four-dimensional de Sitter 
universe emerges purely from quantum fluctuations This is a highly non-trivial 

result keeping in mind that one is dealing with a background independent formulation. 
Other approaches to quantum gravity include string theory, loop quantum gravity and 
causal sets but it is the results on the nature of space-time obtained using the CDT and 
fixed point calculations that particularly concern us here. 

To discuss the nature of a quantum space-time at any distance scale requires a 
quantitative measure which is universal; that is to say it can be defined in any of the 
models we are interested in and is insensitive to cut-off scale physics while conveying 
information about the longer distance structure. The simplest such characterisations are 
various definitions of dimension of which the most familiar is the Hausdorff dimension. 
The Hausdorff dimension d% is defined provided the volume V(R) of a ball of radius R 
takes the form 



if R is large enough. An alternative probe of structure comes from the behaviour of 
unbiased random walks (equivalently diffusion) in whatever ensemble of space-times is 
being considered. The probability P(t) that the walk returns to its starting point after 
time t provides one of the simplest probes of the nature of space time in quantum gravity 
models. The spectral dimension d s is defined if 



For which values of t this should be true is a subtlety with which we will be concerned in 
this paper. For random walks on discretised graphs where the walker is allowed to hop 
from one vertex to a neighbour at each time step, behaviour of the form ([2]) is expected 
at t — y oo i.e. when the walk is much longer than the short distance cut-off scale. On 
the other hand in the continuum the classic picture is that ([2]) describes the behaviour 
as t — > 0. Of course these are not, at least in principle, incompatible as they can be 
related by a rescaling. 

An unexpected result from the numerical simulation of CDTs is that the spectral 
dimension apparently varies from four at large scales to two at small scales [9]. 



V(R) ~ R dh 



(1) 



p(t) ~ r ds/2 . 



(2) 
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Very recently similar results for this phenomenon of dimensional reduction have also 
been observed by other approaches such as the exact renormalisation group JTUJ ITT] . 
Horava-Lifshitz gravity [T2l [13] , and in three-dimensional CDT [H] and some further 
implications are discussed in (15j. Such behaviour is not a priori implausible in 
quantum gravity because there is a dimensionful parameter, namely Newton's constant 
or equivalently the Planck length. Studying the spectral dimension at different 
distance scales raises questions of definition and in the case of numerical simulations, 
discretisation problems. In a numerical simulation the largest available distance scale 
is determined by what will fit in the computer and short distance scales are often not 
much greater than the ultraviolet cut-off, or discretisation scale. Ideally there should 
be a hierarchy in which the long distance scale is much greater than the short distance 
scale which is in turn much greater than the cut-off. In this paper we develop a simple 
model based on previous work on random combs [3]. These are a family of simple 
geometrical models which share some of the properties of the CDT model; instead of 
an ensemble of triangulations we have an ensemble of graphs consisting of an infinite 
spine with teeth of identically independently distributed length hanging off (we define 
these graphs precisely in Section 2). It was shown in [3] that the spectral dimension 
is determined by the probability distribution for the length of the teeth. In this paper 
we show that it is possible to extend the work of [3j by taking a continuum limit thus 
ensuring that the cut-off scale is much shorter than all physical distance scales. We 
find that the spectral dimension is one if we take the physical distance explored by the 
random walk to zero and there exists a number of continuum limits in which the long 
distance spectral dimension differs from its short distance counterpart. As a by-product 
of this work we also extend some of the proofs given in [3] to a wider class of probability 
distributions. 

This paper is organized as follows. In Section 2 we briefly review some known 
results for combs and their spectral dimension and then explain how in principle these 
can be extended to show different spectral dimensions at long and short distance scales. 
In Section 3 we introduce a simple model which we prove does in fact exhibit a spectral 
dimension that is different in the UV and IR. This model forms the basis of all later 
generalisations. In Section 4 we generalise the results of Section 3 to combs in which 
teeth of any length may appear with a probability governed by a power law. In Section 
5 we examine the possibility of intermediate scales in which the spectral dimension 
differs from both its UV and IR values. In Section 6 we analyse the case of a comb in 
which the tooth lengths are controlled by an arbitrary probability distribution and show 
that continuum limits exist in which the short distance spectral dimension is one while 
the long distance spectral dimension can assume values in one-to-one correspondence 
with the positions of the real poles of the Dirichlet series generating function for the 
probability distribution. We then show how these techniques can be used to extend the 
results of [3]. In Section 7 we discuss our results and possible directions for future work. 
Some technical matters are contained in the appendices. 
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Figure 1. A comb. 



2. Combs and Walks 

In this section we review some basic facts about random combs and random walks. As 
much as possible we use the same notation and conventions as [3] and refer to that paper 
for proofs omitted here. 

2.1. Definitions 

We use the definition of a comb given in [3]. Consider the nonnegative integers regarded 
as a graph, which we denote N^, so that n has the neighbours nil except for which 
only has 1 as a neighbour. Furthermore, let Ni be the integers 0,1, ... ,£ regarded as a 
graph so that each integer n G Ni has two neighbours nil except for and £ which only 
have one neighbour, 1 and £—1, respectively. A comb C is an infinite rooted tree-graph 
with a special subgraph S called the spine which is isomorphic to iV^ with the root at 0. 
At each vertex of S, except the root 0, there is attached one of the graphs Ni or A^. We 
adopt the convention that these linear graphs which are glued to the spine are attached 
at their endpoint 0. The linear graphs attached to the spine are called the teeth of the 
comb, see figure [TJ We will find it convenient to say that a vertex on the spine with no 
tooth has a tooth of length 0. We will denote by T n the tooth attached to the vertex 
n on S, and by Ck the comb obtained by removing the links (0, 1), . . . ,{k — 1, k), the 
teeth Ti, . . . ,Tj~ and relabelling the remaining vertices on the spine in the obvious way. 
The number of nearest neighbours of a vertex v will be denoted o~{v). 

It is convenient to give names to some special combs which occur frequently. We 
denote by C = * the full comb in which every vertex on the spine is attached to an 
infinite tooth, and by C = oo the empty comb in which the spine has no teeth (so an 
infinite tooth is itself an example of C = oo). 

Now let C denote the collection of all combs and define a probability measure v 
on C by letting the length of the teeth be identically and independently distributed by 
the measure «. We will refer to the set C equipped with the probability measure v as a 
random comb. Measurable subsets A of C are called events and v(A) is the probability 
of the event A. The measure of the set of combs A with teeth at n\, n 2 , . . . , having 
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lengths £i,£ 2 , ■ ■ ■ ,4 is 

k 

v(A)=n^). (3) 

i=i 

For any z/-integrable function F defined on C we define the expectation value 



(F(C)) = J F(C)du. (4) 
We will often use the shorthand F for (F(C)). 

2.2. Random Walks 

We consider simple random walk on the comb C and count the time t in integer steps. 
At each time step the walker moves from its present location at vertex v to one of the 
neighbours of v chosen with equal probabilities <t(i>) _1 . Unless otherwise stated the 
walker always starts at the root at time t — 0. 

The generating function for the probability pc(t) that the walker is at the root at 
time t, having left it at t — 0, is defined by 

oo 

Qc(x) = J2^-xf 2 Pc(t) (5) 

t=0 

and we denote by Pc( x ) the corresponding generating function for the probability that 
the walker returns to the root for the first time, excluding the trivial walk of length 0. 
Since walks returning to the root can be decomposed into walks returning for the 1st, 
2nd etc time we have 

It is convenient to consider contributions to Pc{x) and Qc{x) from walks which are 
restricted. Let Pq\x) denote the contribution to Pc(x) from walks whose maximal 
distance along the spine from the root is n and define 

n-l 

P { < n \x) = Y,P { c\*) (7) 

k=0 

which is the contribution from all walks which do not reach the point n on the spine. 
Similarly we define 

oo 

P^s) = (8) 

k=n 

Clearly Pc(x) can be recovered from P^ n \x) by setting n — > oo. We define the 
corresponding restricted contributions to Qc{x) in the same way. By decomposing 
walks contributing to Pj^ n \x) into a step to 1, walks returning to 1 without visiting 
the root, and finally a step back to the root it is straightforward to show that 

p c <n \ x ) = ~ — — p «n-i)/ :> ( 9 ) 
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where we have adopted the convention that for the empty tooth, T = 0, 

P $ {x) = 1. (10) 
The relation ([9]) can be used to compute the generating function explicitly for any comb 



with a simple periodic structure and we list some standard results in Appendix A 

There are a number of elementary lemmas which characterise the dependence of 
Pc{x) on the length of the teeth and the spacing between them [3]. We state them here 
in a slightly generalized form which is useful for our subsequent manipulations. 

Lemma 1 The function P^ n \x) is a monotonic increasing function of Pr k {x) and 
Pq^ 1 k \ x ) f or an V n > k > 1. 

Lemma 2 Pj^ n \x) is a decreasing function of the length, 1^, of the tooth for any 
n > k > 1. 

Lemma 3 Let C be the comb obtained from C by swapping the teeth Tk and T^+i, 



k < n — 1. Then Pq\x) > P^ n \x) if and only if PT k (x) > PT k+1 i 



x 



The proofs use ([9]) and follow those given in [3] for the case n = oo. 

An important corollary, valid for any comb, of these lemmas is that 

< Q c (x) < x~^, (11) 

which we will refer to as the trivial upper and lower bounds on Qc{%)- The result follows 
from Lemma [2] with n = oo, which gives 

P*{x) < P c {x) < Poo{x), (12) 



and the explicit expressions for P*(x) and Poo(x) given in Appendix A 



2. 3. Two point functions 

Two point correlation functions on the comb correspond to the probability of a walk 
beginning at the root being at a particular vertex on the spine at time t. In particular, 
let pc(t;n) denote the probability that a random walk that starts at the root at time 
zero is at the vertex n on the spine at time t having not visited the root in the intervening 
period. We will refer to the generating function for these probabilities as the two point 
function, Gc(x;n), and define it by 

oo 

G c (x;n) = Y,(~L-x) t/2 Pc(t;n). (13) 
t=i 

Gc(x; n) may be expressed as 

n-l 

Gdx; n) = a(n)(l - x)^ 2 J] P Ck {x) (14) 

k=0 

which may be used in conjunction with Lemma [2] to obtain the bounds, 
G*(x;n) < G c (x;n) < G^^n) 

3 ~ a(n) 2 1 ' 
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Now let rc{t;n) denote the probability that a random walk that starts at the root at 
time zero is at the vertex n on the spine for the first time at time t having not visited 
the root in the intervening time. We define the modified two point function, G^(x;n), 
by, 

oo 

G c (x- 1 n) = Y,(l-x) t/2 rc(t ] n) (16) 
t=i 

and note the following lemmas; 

Lemma 4 The contribution Pq (x) to Pc{x) from walks whose maximal distance 
from the root is N or greater satisfies 

P [ c N ' l \x) < 3x-^ 2 G° c (x; AO 2 . (17) 
The proof is given in section 2.4 of [3J. 
Lemma 5 The modified two point function satisfies 

G^(x;n) <G§\x;n) <Gg?(s;n). (18) 
To prove this note that 



{X, lb) — 

and use Lemma [2] 



^ > k=0 



2.4- Spectral dimension and the continuum limit 

The spectral dimension of random combs was studied in [3] . In this work the probability 
distributions /jl(£) for the length £ of a tooth were chosen to be fixed sets of numbers 
so the teeth at adjacent sites on the spine are not only independent but generally show 
large fluctuations relative to each other. In these circumstances the spectral dimension, 
d s , describes the large t dependence of the return probability pc(t) and is given by 

d s = -2 lim lQg ^ (t)) (20) 

t^oo log t 

if the limit exists. In fact it is much more convenient to deal with generating functions 
and by a Tauberian theorem [17] we expect that if (|20|) holds then, as x — > 0, 

Q c (x) ~ x- 1+d °/\ (21) 

where by f(x) ~ g(x) we mean that 

cg(x) < f(x) < c'g(x), < x < x , (22) 

where c, d and xq are positive constants. The property ( I2T1) was adopted in [3] as the 
definition of spectral dimension, assuming it exists. Heuristically the spectral dimension 
characterizes certain aspects of the long distance structure of a graph as observed by 
a walker who goes on a very long walk and hence probes that structure. The spectral 
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dimension of an ensemble average is defined in the same way, simply replacing pc{t) 
and Qc(x) by their respective expectation values. 

In this paper we study the possibility of different spectral dimension on different 
distance scales. To do this we have to generalize our definition from (120]) or (I2TT) and 
introduce at least one characteristic distance scale L f into the probabilities fi(£) 
which determine the structure of the comb. We then assign the value a to the distance 
between adjacent vertices in the graph and take the limit a — > and L — > oo in such 
a way that the scaled combs have a finite characteristic distance scale; we will refer to 
this limit as the 'continuum' limit and quantities which exist in this limit as continuum 
quantities. Walks much longer than L will probe different structure from walks much 
shorter than L but nonetheless both can be very long in units of the underlying cut-off 
scale a. 

In the following sections we will denote dependence of a function on a number of 
variables Lj, i = 1, . . . ,N, by Lj passed as one of the function arguments. Given a 
random comb ensemble specified by /i(£; Li) and the corresponding Q(x; Li) we define 

K) = lim <AQ(a£;a- A *Af<), (23) 

a— >-0 

where the scaling dimensions A M and Aj are chosen to ensure a non-trivial limit and 
the combinations £Aj are dimensionless. Q can be used to define the spectral dimension 
at short and long distances. 

In the following discussion we assume for simplicity that there is just one scale L 
and that the spectral dimension in the sense of (l20i) exists for (pc{t)) which implies that 
there exists a constant t such that (pc(t + I)) < (pc(t)) ,t > t . Note that 

T oo 

J2(Pc(t))(l-xf 2 =Q( X] L)- (Pc(t)) (1 - xf 2 

t=0 t=T+l 

vt/2 



Q(x; L) - (1 - xY T+1 » 2 J2 (Pc(t + T + 1)) (1 - xf 



t=o 

oo 

> Q(x; L) - (1 - xY T+1 ^ 2 <Pc(t)) (! - ^ ■ (24) 

i=0 

Now choose 



T 



a- 1 



£iog(i + s) 



1 (25) 
and set x = a£ and L = a~ A \ A in ( |2^1) to get 

T 

a A »Q{a£; a~ A \ A ) (1 - exp (-£ A)) < a A » ^ (p c {t)) (1 - ^)* /2 < a" A A A !j26) 

Provided that the limit in fl23l exists we see that the behaviour of A) as £ — > oo 
characterizes the properties of walks of continuum time duration less than 

^oo^log(l + 
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and we define the spectral dimension d° at short distances by 

iog(Q(& 

hm — - 

C^oo log 4 



d° = 2[l+lim 1 ° g( 1 ^; A)) l, (2f 



provided this limit exists. 

We can define the spectral dimension at long distances in a similar way. First note 
that by §TB) 

f 1 \ 1,2 T f 1 \ t/2 f 1 \ T T 

so that 

oo 

Q(x; L) - VT(1 - -Y T < J2 (Po(t)) (1 " *f 2 < L). (30) 

t=T+l 

This time letting T = [a -1 ^ -1 log(l + £A)J — 1 we get 

oo 

cAQ«; a~ A A A ) - ey^" 1 log(l + JX) < a A " ]T (p c {t)) (1 - ^af 2 < <AQ«; a" A A A ). 

t=T+l 

(31) 

Provided that the limit in (1231) exists and that <3(£; A) diverges as £ — >■ we see that 
its behaviour characterizes the properties of walks of continuum time duration greater 
than lim^o £ _1 log(l + £A) = A . We then define the spectral dimension df at long 
distances to be 

dr^fi + 1 . m ' og ff> ) V P2) 

^ log£ J 

provided this limit exists. 

It is by no means obvious that there are graph ensembles for which the limits (I2"3"j) 
followed by f )28|) and f l32j) exist. However in the rest of this paper we will show for 
comb ensembles of increasing generality that this is indeed the case. Clearly at the very 
least any such ensemble must have a characteristic distance scale A that survives the 
continuum limit otherwise such behaviour is impossible. In all the examples given in 
this paper it turns out that the exponent A M = ~. 

3. A simple comb 

We now introduce a random comb whose spectral dimension differs on long and short 
length scales and thus illustrates that the behaviour described in section 12^1 can actually 
occur. This comb is defined by the measure, 



1-1 i = 0, 

7, £ = oo, ( 33 ) 

0, otherwise. 
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This random comb has infinite teeth and they occur with an average separation of 
L. Intuitively we would expect that if a random walker did not move further than a 
distance of order L from its starting position it would not see the teeth and therefore 
would measure a spectral dimension of one. If however it were allowed to explore the 
entire comb it would see something roughly equivalent to a full comb and so feel a much 
larger spectral dimension. To prove this intuition correct we proceed by computing 
upper and lower bounds for Q which are uniform in L and for < x < Xq, where the 
constant Xo is equal to one unless otherwise stated, and then take the continuum limit 
to obtain bounds for Q . 

With complete generality we may obtain a lower bound on Q(x) by use of Jensen's 
inequality which takes the form, 

Lemma 6 Let Pt{x] Li) be the mean first return probability generating function of the 
teeth of the comb defined by \x{l\L^), then 

Q(x;Li) > (1 + x - P T (x; Li))-? . (34) 

The proof is given in [3]. For the comb ( 1331) we have 



P T (x; L) = 1 - 1(1 - P^x)) = 1 - Y (35) 



which implies 



Q{x;L)>i^ + x) . (36) 
Letting x = at; and L = a~z\^ gives 

_ i 

Q(f; A) = lim alQ«; arh?) > (v~^ + l) \ (37) 

where we have introduced the dimensionless variable v = £A. 

To find an upper bound on Q(x; L) we follow [3] and use Lemmas [fl [2] and [3] to 
compare a typical comb in the ensemble with the comb consisting of a finite number of 
infinite teeth at regular intervals. First we define the event 

A(D, k) = {C : Di < D : i = 0, k}. (38) 

where is the distance between the % and i + 1 teeth and then write, 



Q(x;L)= J Q c (x;L)du 



Q c (x;L)dis+ / Q c (x;L)dis. (39) 

C/A{D,k) JA(D,k) 

Since the Di are independently distributed 

u(A(D,k)) = (l-(l-l/L) D ) k . (40) 
Consider a comb C G A(D, k); then by Lemmas [fl [2]and[3l 

P c (x;L)<P c ,(x), (41) 
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where C is the comb obtained by removing all teeth beyond the k tooth and moving 
the remaining teeth so that the spacing between each is D. Now we can write 

P c ,(x) = P^ Dk \x) + p£ Dk - 1] (x). (42) 

Since the walks contributing to P^ Dk \x) do not go beyond the last tooth we have 

P£ Dk) (x) < P, D (x), (43) 

where *D denotes the comb consisting of infinite teeth regularly spaced and separated 
by a distance D. Using ()43l) . Lemmas H] and [5] we have, 

P c (x; L) < P* D (x) + 3x^G^\x; Dkf (44) 

uniformly in A. P*d{x) and G^(x; n) are given in Appendix A. Now set D = \_D\ and 
k = \k], where, 

D = 2L\ logxL 2 |, k = (xL 2 y 1/2 . (45) 

Since G^(x;n) is manifestly a monotonic decreasing function of n and P*d{x) an 
increasing function of D, 

Q(x; L) < x- l '\l - (1 - (1 - 1/Lf -Y+ 1 ) + Q v (x)(l - (1 - l/Lff (46) 
where we have used ffTTj) and 

Qu(x) = fl - P^(x) - 3x- l 2G<g(x; (D - l)~k) 2 } ^ . (47) 



Taking the continuum limit of (T4"6"j) and using the results of Appendix A then gives 

g(£;A)<r 1/2 i^A), (48) 

where 

1 + o(v _1 ), v — >• oo, 

F(v)={ , ' _ , (49) 

y I logf 2 | + o(t> 2 ), v — > 0. 



It follows from (128lh 032]), (I37j) and (148]) that 

<£ = 1, d? = | (50) 

4. Combs with Power Law Measures 

We now consider slightly more general combs in which the measure on the teeth is a 
power law of the form, 

fl-y, 1 = 0, 

^■L)= I (51) 

[ -c a r a , i > 0, 

where C a is a normalisation constant and as before L plays the role of a distance scale. 
We consider laws in the range 2 > a > 1 as it is known that for a > 2 the comb has 
spectral dimension d s = 1 in the sense of (12T!) [3] and therefore it is not possible to get 
a spectral dimension deviating from 1 on any scale. 
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To compute a lower bound on the return probability generating function for the 
above distribution we apply Lemma [6] and reduce the problem to computing an upper 
bound on 1 — Pt(x). The first return generating function Pi{x) for a tooth of length I 
is recorded in (1A.5I) ; bounding tanh(w) above by the function f(u) = u for u < 1 and 
f(u) — 1 for u > 1 give^~" 



l-P T (x; < 



L ^ + ^ u 



p ^T-r / M \ 

/ moc{x) I _ — ^ I 

< a/x - rriooia^y/x I \ //(£;Lj) cfat. (52) 



To obtain the second inequality we have applied the Abel summation formula. We 
therefore have, 



m 00 (x) 

x{u; Li)du (53) 



Lemma 7 For a random comb defined by the measure \x, 

i 

1 - P T {x] Li) <\fx- rrioo^y/x / 

Jo 

where the cumulative probability function x(w; Li) is defined by x(«; Li) = 2^[=o M^j ^i)- 

We will see shortly that all behaviour of the spectral dimension of the continuum comb 
is encoded in the asymptotic expansion of x( u 'i Li) as u goes to infinity 0. In the present 
case x{ u \ L) is trivially related to the partial sum of the Riemann ^-function whose 
leading asymptotic behaviour is well known and we find 

x ( u -L) = l-^^—- + 5(u), (54) 
L a — 1 

where 

\5(u)\ < yu- a , u>2. (55) 

LJ 

It follows that for x < xq, where moo(a;o) = |, 

1 - P T (x; L) < moo (x)V^ (^j- moo (x) a - 2 + b j- moo {x) a - 1 + , (56) 

with 61,2,3 being constants depending only on a and 61 > 0. Choosing L = a~ A ' X A ' with 
A' = 1 — a/2 yields a lower bound on the continuum return generating function, 

Q^x)>r 1/2 (i + hm- {i - a/2) y 1/2 . (57) 

To obtain a comparable upper bound we need 

I In this particular case we could in fact compute 1 — Pt{%) exactly by the Abel summation formula. 
However the bound we use is good enough to give the desired result with the advantage that the 
calculation can be done with elementary functions. 

§ In general it is not obvious that this asymptotic expansion exists due to the discontinuous nature of 
X- We will address this issue later when we consider generic measures. 
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Lemma 8 For any random comb and positive integers H , D and k, the return 
probability generating function is bounded above by 

Q(x; < x- x l\\ - (1 - (1 - p) D ) k ) + Qu(x)(l - (1 - p) D )\ (58) 



where 



e=H+i 

Qu{x) = [l - PhM*) ~ 3aT3G2?(a;; Dk) 2 j ' , (59) 

and Ph,*d is the first return probability generating function for the comb with teeth of 
length H + 1 equally spaced at intervals of D. 

The proof is a slight modification of the upper bound argument used in Section [31 First 
define a long tooth to be one whose length is greater than H; then the probability that 
a tooth at a particular vertex is long is 

oo 

P= Yl (60) 

i=H+l 

Define the event 

A(D,k) = {C : A < D : i = 0, ...,*} (61) 
where now Di is the distance between the i and i + 1 long teeth so that 

Q(x;Li)= j Q c (x;Li)du 



Qc(x; U)dv + / Qc(x;Li)du. (62) 

C/A(D,k) JA(D,k) 

Since the Dj are independently distributed 

K^(AA:)) = (l-(l-j9) D ) fc . (63) 

Now use Lemmas [2] and [3] in turn to note that for 

PczA{D,k){x,L)<P c ,{x,L) (64) 

where C is the comb in which all teeth but the first k long teeth have been removed and 
the remaining long teeth have been arranged so that they have length H and a constant 
inter-tooth distance D. By the same arguments as we used in Section 3 to get (1541 we 
obtain the bound 

PczA{D,k)(x,L) < P H ^ D + 3x~ 1/2 G { °\x,Dk) 2 . (65) 

Lemma [8] then follows from (jlip . (1631) and (j63j) . We now specialise to the power law 
measure ( 15"T1) and set if = \_H\ , D = [-DJ and k = \k], where 

H = x- 1 / 2 

D = (A' + l)^!^'- 1 / 2 ^ logxL 1 ^'! (66) 
= (xL 1 ^')-* 
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Using Lemma El the scaling expressions for Ph,*d and given in ( 1A.7I) . and 
taking the continuum limit, gives, after a substantial amount of algebra, 

Q&X)<C 1/2 F(Z\), (67) 

where 

fl + Oftr 1 ), v ->■ oo, 

= <^ , V A , (68) 

\c^ 2 - Q / 4 v^loi^+0(t; A ), u->0. 

The main result of this section is 
Theorem 9 TTie comb with the power law measure (51\) for the tooth length has 

d° s = 1, rff = 2 - |. (69) 
The result follows immediately from Q2g), (I52|l . (]5Tj) and (]oTj). 

5. Multiple Scales 

Given the results for the power law distribution it is natural to investigate the behaviour 
for a random comb that has a hierarchy of length scales. The easiest way to achieve 
such a comb is through a double power law distribution, 

( -L 2 \ £ = 0, 

n(£- LA = < 1 1 (70) 

We may assume without loss of generality that the length scales Li scale in the continuum 
limit to lengths Aj such that Ai < A 2 and that 1 < < 2. 

Following the procedure of previous sections a lower bound on Q(£; A«) is obtained 
by using Lemma [6] and noting that x{ x > Li) for this comb is essentially the sum of the 
cumulative probability functions for each power law. This gives 

1 - Pt(x; U) < moo (x)V^J2 (¥™<oc(xp- 2 + j^m^xp- 1 + h f\ . 

i=i ^ 1 1 1 ' 

(71) 

Choosing Li to scale like Li = a _A *Af% where = 1 — an/2, gives a bound on the 
continuum return generating function of, 

T 1/2 (co + c^XJ-^M + c 2 ^X 2 )^-^y 1/2 < Q(£). (72) 

An upper-bound on Q(£) is obtained by application of Lemma [8] in which we set 
H = [D\ , D = [D\ and jfe = \k], where 

l) = /Jx-^G^L^SxL^)- 1 ! logxL^' 1 ! (73) 
k =G(xL\ /A '\xLl /A ' 2 ) 
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and for convenience we have introduced the function 



r*t \ ^i -A' C 2 —A' 2 

G(v 1 ,v 2 ) = -v x l + -v 2 2 . (74) 

«i — 1 a;2 — 1 



Using Lemma |8] and the scaling expressions in Appendix A then gives 
Q(x;A,)<r 1/2 [l-(l-"^ G 

(i-<7 



(1 - v- sl T + 



(75) 



3cosech 2 (|logvf|(l - l/G)) - 7 + ^/ 7 2 + 1 + 2 7 cothfl log |/G)- 

where u$ = £Aj, s = sgn(log-Ui), 7 = tanh(l) and we have suppressed the arguments of 
G(vi,v 2 ) in order to maintain readability. 

We can now examine (172]) and (1751) to see what they tells us about the behaviour 
of Aj) on various length scales. 

• When £ ;§> A]" 1 both upper and lower bounds of Aj) are dominated by the 
behaviour so taking the £ — >■ 00 limit leads to = 1 as in the previous sections. 

• If ct\ < a 2 then when £ <C A^f 1 both upper and lower bounds of <3(£;Aj) are 
dominated by the £ -ai / 4 behaviour so taking the £ — >■ limit leads to df = 2—a\/2. 
There is no regime in which a 2 controls the behaviour. 

• If a 2 < ct\ then when £ <C A -1 where 

both upper and lower bounds of Q(£; Aj) are dominated by the £ _Q2 / 4 behaviour so 
taking the £ — > limit leads to df = 2 — 02/2. However there is an intermediate 
regime A -1 <C £ <C A^f 1 where the £~ ai / 4 behaviour dominates and Aj) lies in 
the envelope given by 

Cl r Ql/4 < A,) < c 2 r Ql/4 v^g^ ; (77) 

where the upper and lower bounds will have corrections suppressed by powers of 
£\ 2 and the upper bound will also have corrections of order £^. Both A2 and (5 may 
be chosen to make the corrections arbitrarily small in this scale range. The system 
therefore appears to have spectral dimension 5 s = 2 — a±/2 in this regime. This 
is a fairly weak statement because Q(£) could in principle exhibit a wide variety 
of behaviours between its upper and lower bounds; this region is just a part of the 
crossover regime from d s to 4%°. However, as we are free to chose A2 to be as large 
as we like compared to Ax, this regime can exist over a scale range of arbitrarily 
large size. We therefore can force the leading behaviour of Q(£; Aj) in this range to 
be as close to a power law with exponent 5 s = 2 — cti/2 as we like. This is what 
might be observed, for example, in a numerical simulation; if the difference between 
the scales Ai and A2 is large then there will be a substantial range of walk lengths in 
which the data will indicate a spectral dimension of 5s- We will refer to a spectral 
dimension that appears in this weaker way as an apparent spectral dimension and 
denote it by 5s rather than ds- 
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6. Generic Distributions 

So far we have considered combs in which the distribution of tooth lengths has been 
governed by power laws or double power laws. In this section we extend the results of 
the previous sections to the case were the form of the tooth length distribution is left 
arbitrary. The most general situation is that the measure on the combs is a continuous 
function of some parameters Wj. The continuum limit of such a comb is obtained in 
the usual way but with parameters wi scaling in a non-trivial way; Wi = w Ci + a dt Ui. 
Given a random comb with such a measure we would like to know how many distinct 
continuum limits exist and for each compute how the spectral dimension depends on 
the length scale. 

The approach we adopt here closely mimics the arguments of the preceding sections, 
indeed the main complication is technical. As we have seen the properties of the 
continuum comb are controlled by the asymptotic expansion of x{ u ) as u goes to infinity. 
The main difference in the generic case is that we may arrange matters so that the 
scaling dimensions of the coefficients in the asymptotic expansion are such that sub- 
leading terms appear in the continuum. For the generic case we obviously have no 
way of knowing the full asymptotic expansion of x{ u 'i w i)- However, we will see that 
for a large class of measures, the form of the asymptotic expansion is encoded in the 
asymptotic expansion of a particular generating function for \i{£\ w^). 

Our first task is to introduce this generating function and relate it to the asymptotic 
expansion of x{ u 'i w i)- To this end we introduce the notion of a smoothed sum [IB"] . 

oo oo 

X±(u;wi) = ^2n(e;wi)r] ± (e/u) = //(O;^) + ^2fi(£;Wi)r] ± (£/u) = /x(0) + x± \w,Wi), 
e=o 1=1 

(78) 



where rj± is the smooth cut-off function introduced in Appendix B and u controls where 
the cut-off occurs. Such smoothed sums are related to x{ u \ w i) by, 

X-(u; Wi) < x(u] Wi) < x+(u; w^. (79) 

The reason for introducing the smoothed sums is that we may use powerful techniques 
from complex analysis to compute their asymptotic expansion (see eg [17]). The 
generating function to which the asymptotic expansion of the smoothed sums is related 
is the Dirichlet series generating function of fi(£;Wi), 

w i) = 2^ — T, — • ( 80 ) 

We now introduce a number of results and notations, 

1. For a strip in the complex plane S(6, a) = {z : a < ~Re[z] < b} where b > a, we say 
Vfj,(s;Wi) has slow growth in E(6, a) if for all s G S we have "D M (s;Wj) ~ 0(s r ) for 
some r > as Im[s] — >■ ±oo. We say V^s^Wi) has weak slow growth if the above 
property only holds for a countable number of horizontal regions across the strip. 
See figure [21 
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Figure 2. An illustration of the strip S in which T>^(s; wi) satisfies property (f ) given 
given in the text. The poles of T>^{s\Wi) are denoted by crosses and the horizontal 
strips in which the growth condition holds are indicated by the dark grey regions. 



2. Define S± to be the set containing the triples (— o + it, —k, —r±) such that the 
Laurent expansion of T>^(s; Wi)Ai[^/± t ](s + l)/s about the point —a + it contains 
the term r±/{s + cr — ir) k where k > and —a + it ^ 0. Here M. denotes the 



Mellin transformation and \I/± e is introduced in Appendix B We will often find it 
useful to refer to the elements of S± using an index j and denote the jth element 
of S± by (— Oj + iTj, —kj, — Tj t ±). Note that since A^[^± e ](s + 1) is analytic in E 
then the positions of the poles are determined only by T>^(s\ wi). 

3. Define the indexing sets Jr and Jc, such that if j G Jr then Tj = whereas if 
j G Jc then Tj ^ and for both —Gj + iTj G S i.e. they index the poles in E which 
lie on the real line and off the real line respectively. 

4. For a Dirichlet series with positive coefficients the abscissa of absolute and 
conditional convergence coincide. Furthermore, since YaIo w i) = 1 ^ ne abscissa 
of convergence is less than zero. 

5. Landau's theorem: For a Dirichlet series with positive coefficients there exists a 
pole at the abscissa of convergence. A corollary of this is that the coefficient of 
the most singular term in the Laurent expansion about the abscissa of convergence 
is positive and furthermore it is the right-most pole of the Dirichlet series in the 
complex plane. 

If T>^ is of slow growth in E(a, b) and no pole occurs to the right of E(a, b), then 
by using the expression (1C.10I) the asymptotic expansion of X±{ u 'i w i) 

X±(u; = /i(0; w { ) + ^— <fc w s X> M (s; Wi) M ^ ±€ S + 1 ds - R(u; Wi) 

Ztii J c s 
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1 - > J r ji± W \, ' - (81) 



^ E { r i,±{ w i) uiTi + r h± \w^) ( ^ZW ^ - R{u- Wi) 



where C is the rectangular contour introduced in Appendix (jj j], r^± are the coefficients 



of the Laurent expansions of V^s; Wi)M.[^ ± € ](s+l) / s, and R(u) is a remainder function 
which satisfies R(u) ~ m _7V , where A" > aj, as w goes to infinity. Note that the 
difference between the coefficients rj t+ and rj_ is of order e, the parameter introduced 
in 



Appendix B , which can be taken to be arbitrarily small and so r Jrf and r 3 - are for 
all purposes equal. 

We now are in a position to prove the main result of this section. Note that in the 
following we require the continuum comb is such that A M = 1/2 and so has spectral 
dimension one on the smallest scales. It seems very likely that this restriction could be 
lifted but we will not pursue such a generalisation here. 

Theorem 10 For a comb in which the teeth are distributed according to the measure fi 
then if Up has slow growth in the strip S = {z : — 1 < Re[z] < 0} and has no poles on 
its line of convergence besides at the abscissa then continuum limits of the comb exist 
with d® — 1 and df taking values in the set {^-p- : < Oj < l,j G Jr}- 

Proof The proof proceeds in much the same way as the proofs in the previous sections; 
we first derive upper and lower bounds on Qc and then use these bounds to deduce the 
behaviour of Qc in different scale ranges. 

We begin with the lower bound which by Jensen's inequality amounts to finding 
an upper bound on 1 — Pp. Note that implies that the scaling dimension of 1 — Pt 
must be greater than or equal to one in order for A M = 1/2 and only contributes to the 
continuum limit if it has scaling dimension one. From Lemma [7] we have, 



— / "loo (a:) 

P T (x; Wi) < m 00 (x)Vx / (1 - X-W wi))du 

Jo 

i 

/ m M (i| 

< y^xm^^C^; Wi) + m 00 (x)\rx / (1 — 

Jk 



X~(u;Wi))du, (82) 



where C(K) = f^(l — X-{ u \ Wi))du is a constant independent of x. It is important to 
recall that since < x±{ u '> w i) — 1 f° r an values of u and W\, . . . , wm, in particular when 
Wi, . . . ,wm assume their critical values, that the scaling dimensions of the coefficients 
r j,± appearing in ( )8T|) must be positive as otherwise x±( u i w i) would diverge if we were to 
set wi, . . . , wm to their critical values. Upon performing the integration in (1521 we will 
find that a given rj ± now is the coefficient for a number of terms of increasing scaling 
dimension. Since we require the scaling dimension of 1 — Pt to be greater than or equal 
to one, only the term with smallest scaling dimension can appear in the continuum 

|| Since we have assumed there are no poles to the right of S we may choose c = b, where c is the 
constant appearing in the definition of the contour. 
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limit and we will drop any term that does not appear in the continuum limit. If we now 
substitute flgTJ int o flH2D we find 

1 - P T {x] Wi) < y/xm OD (x)C(K; wi) + ~ X 




(83) 

where the remainder term R disappears since it cannot appear in the scaling limit. We 
now suppose we may choose the critical values and scaling dimensions of the parameters 
Wi such that Tj ± has a scaling form that can be written as, 

where 9j, 9j and \j t ± are constants and we have included the factor of | for later 
convenience. From (|83|) one can see that since we require A M = 1/2 the only terms which 
appear in the continuum limit are those for which 9j = (1 — o~j)/2 and 9j = — {kj — 1) 
and so it is useful to define a restricted indexing set 

Jr = {] eJ R :6 j = (l-a j )/2 and 9j = -(kj - 1)} (85) 

and an equivalent one Jq for Jq. The continuum limit is thus, 

1 - P T{X ; Wi ) < aCi + £ Ji—J^^i KA*-)-^ + (86) 

E (^T)I ( Re ^] cos (f MO) + Im[0,]sin(|log(aO)^(e|A,,-|)-^ i , 

where 0j = (|Aj i _|/Aj i _)^ 1_f7: '^ 2 (l — Oj +irj)^ 1 . The lower bound on Q c is then obtained 
from Lemma [6j 

Before deriving an upper bound on Qc, we must analyse the consequences of any 
of the oscillatory terms in the second sum appearing in the continuum limit (i.e. if Jc 
is non-empty). In the case of the double power law measure we saw that one could 
obtain intermediate behaviour in which an effective spectral dimension was measured 
that differed from the UV and IR spectral dimensions. A similar phenomenon occurs in 
the generic case when the various length scales of the continuum limit are well separated. 
If we scale the coefficients of the oscillatory terms such that these terms appear in the 
continuum limit then we are lead to an inconsistency by the following argument, 

1. Let the term associated with an oscillatory term have index j . Consider choosing 
the scaling form of x and rj_ such that £|Aj 0( _| 1 but £|Aj,-| 3> 1 if 
|Aj_| > |A jo _|. 

2. If we were to take the scaling limit in such a scenario then the term that dominates 
the size of 1 — Pt is the one associated with the length scale |Aj 0) _|. This term will 
be oscillating around a mean of zero and hence must be going negative infinitely 
often as we approach the continuum. 
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3. Since Pt is a probability generating function it cannot be negative, hence showing 
we have an unphysical limit. 

The cause of this behaviour is that we declared by fiat that the parameters Wi must be 
scaled to give (ISlj) ; however the parameters must also satisfy the constraints Wi) > 
and A*(^> w i) = 1 f° r a ^ va l ues of the parameters and we have not ensured these 
constraints are compatible with It is sufficient for our purposes to understand that 
the oscillatory terms prevent the continuum limit being taken for certain walk lengths 
so they are certainly unphysical; we therefore must scale them so that they disappear 
in the continuum. We therefore have as a lower bound on Q c , 

-1/2 



<?(£; <*) > r 1/2 i + s AC{K) ,oC + £ {k . _ 



-v 2 



°3, 



(87) 

We now consider the upper bound. Without loss of generality we may arrange that 
the indexing set Jr has the property that if j\ < j 2 then Xj lt ± < Xj 2i ±. Furthermore 
define 

H(x) = [H(x)\ = [x- l/ \ 

D(x) = [D(x)\ = - X (^)))- 1 |log(r /i+ x^(-ilogx)^- 1 )|], (88) 

I — o~j I 

k(x) = [k(x)} = [x- 1 / 2 (l~x(H(x)))], 

where I = inf J^, i.e. Xj t ± is the smallest length scale with a spectral dimension differing 
from one. Some modified versions of the above quantities will also be needed, 

D ± (x) = [D ± (x)\ = [-^(1 -x ± (^)))- 1 |log(r / , + x^(-^logx)^- 1 )|], (89) 

1 — a j I 

k ± (x) = [k ± (x)} = [x- 1 / 2 (l- X ±(H(x)))]. 

It is clear that D_(x) < D(x) < D + (x) and k + (x) < k(x) < k_(x) which together with 
Lemma [8] allows us to conclude that the continuum limit is, 



Q(x;ui) <r V2 



1 _ e G_log(l-«7 3 ' 3 ) + 



e (G+-l)log(l-«7 s ' 3 ) 



3cosech'(| log u7 1(1 - 1/G+)) - 7 + a/7 2 + 1 + 27coth(| logvf |/G+) J 



(90) 



where we have defined G±(£, Ai, . . .) = ^2j e j R Vj t ± 2 /(kj — 1)! and arranged that no 
oscillatory terms appear. 

We are now in a position to prove Theorem [10] by analysing the behaviour of 



( I87|) and (190]) for various walk lengths. We first note that both the upper and lower 
bounds are controlled by the relative sizes of the quantities i/ . 1+<Tj / 2 . i n particular the 
largest of these quantities, T = sup{u 1+ct j/ 2 : j g J R y^ w m determine the behaviour 
in a particular scale range; if we suppose T = v T 1+aj/2 for some j then the leading 
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contribution to both the upper and lower bound will be proportional to £ ( 1+ °j)/ 4 . 
Using this we find the following behaviour, 

• On very short scales corresponding to £ ^> Xj 1 for all j then the lower bound of 
<5(£) is dominated by the behaviour, which together with the trivial upper 
bound means that taking the £ — > oo limit leads to d° = 1 as in the previous 
sections. 

• On very long scales corresponding to £ Xj 1 for all j then there will exist a length 

scale A such that for < £ < A -1 , T = v- + ^ where a- < o~j for all j. Explicitly 
A will be given by, 

A = supjAj- A, :je J R \. (91) 

If we take the limit £ — >■ we obtain = (3 — erj)/2. 

We see that the spectral dimension increases monotonically on successive length 
scales and the spectral dimension has measured values given by d s — (3 — cr-vn)/2 
therefore proving Theorem [TOl 

It is interesting to note that a constraint on the crossover behaviour exists for generic 
measures much as it did for the case of the double power law measure. In particular on 
intermediate scales there exists J E Jr such that £ <C Xj 1 for j < J and £ ^> Xj 1 for 
j > J. Hence Vj 1+a ^ 2 <c 1 if j > J and so T = sup{^- 1+a ^ 2 ■ j g j R ^j < j}. This 

means there will exist a length scale Aj such that for Aj+i -1 < £ < Aj _1 , T = Vj,l J(J) ^ 

where we have defined j( J) implicitly by < o~j for all j < J. Of course this does not 
ensure that Aj + i _1 < Aj -1 , indeed if the length scales Aj are not sufficiently separated 
this may not be true and we would not have a scale range in which this term dominated. 
The expression for Aj is, 

Aj = supjA^ A, (J) J : j G J B , j < J} (92) 

and so we may choose Aj + i independent of Aj thereby allowing the scale range over 
which this behaviour exists to be arbitrarily large. This would result in an apparent 
spectral dimension of S s = (3 — c--/n)/2. 

Finally, an interesting application of the techniques used to prove Theorem [10] is 
that they allow one to analyse a wider class of combs than the class for which the results 
of [3] are valid. In particular it was proven in [3] that for a random comb the spectral 
dimension is d s = (3 — 7q)/2 where 70 = sup{7 > : J 7 < 00} and J 7 = J2^o t^i ^ ■ 
This was proved subject to the assumption that there exists d > such that, 
00 

E (93) 
as x goes to zero. 

Given the results in this section we see that —70 may be interpreted as the abscissa 
of convergence for the Dirichlet series generating function. Furthermore, it is clear from 
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our results that there are distributions where the assumption (1931) does not hold and 
that we may use the techniques we have developed to analyse these cases. Recalling 
that for a random comb we may compute the spectral dimension using the relation ( 12T1) 
then we must perform a similar analysis to that done for the continuum comb but now 
only scaling x to zero. 

Due to Landau's theorem there always exists a pole at the abscissa. Suppose it is 
of order k and consider 1 — x±( u )i 

-logy) fe -' 



i - x±(y) 



r ,± + ^ Re t r i.±] cos ( r i lo S?/) + Im [ r i,±] sin ( r i lo SV) 



(k 



-V 



10 



tt{y; {r jt ±}) 



-log?/) 



fc-i 



,7() 



(94) 



(k-l)\ • 

where j runs over the poles on the line of convergence and fl(y; {?"j,±}) is implicitly 
defined in the second line. We have not included any poles of order less than k since 
these will not be leading order as y goes to zero and we have assumed the there are no 
poles of order greater than k as then the above quantity would go negative due to the 
oscillating terms. Applying the same argument as we did for the continuum comb we 
obtain the lower bound, 

c'{k-l)\ i 1/2 ^ ^ (1 - fe)/2 



Q{x) > 



_ft(x 1 /2 ;rj ._(l_ 7o + ^.)-i) 



-| logx\ 



X 



70 + 1 
4 



(95) 



where c' is a constant. To obtain an upper bound we apply a very similar argument as 
used for the continuum comb; the only difference being that we choose, 

H =[H} = /3aT 1/2 |logo;|, 

1-1 logx|, 



D 

k 

D± 
k± 



[D] = 0(1 - X {H)) 
[k}=x- 1 / 2 (l- X (H)), 
[D±]=P(l-X±(H)y 1 \ 
[k±}=x^ 2 (l- X ±(H)) 



log a; | 



and compute the behaviour as x goes to zero of G^(x, (D — l)(k 
size of the set A in ( )6T|) . The result is, 

c"(k - 1) 



1)), ^ 



(96) 



and the 



Q(x) < 



.1/2 



1 1/2 f, 


- x l/2 - 


1 1 log 


log x@ 



(l-fc)/2 



logx 



/3|70/2 x -^± 



(97) 



• |loga;^|' '3>+J 

where c" is a constant. We see that we reproduce the result of [3] when k = 1 and there 
are no poles on the line of convergence. If k ^ 1, the spectral dimension is the same 
as the k = 1 case as only logarithmic corrections are introduced. If we allow poles to 
appear on the line of convergence then Q(x; r 3 - ±) will have an oscillating x dependence 
which, if the bounds above are tight enough, would also imply that the functional form 
of Q is not a power law and so the spectral dimension does not exist. Whether the 
above bounds are tight enough to make this conclusion we leave for future work. 
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7. Conclusions 

We have demonstrated that there exist models in which a scale dependent spectral 
dimension can be shown to exist analytically. That we could do this was important 
as up to now all the evidence for scale dependent spectral dimension in CDT has 
been numerical in nature and therefore open to the criticism that it might be due to 
discretisation effects or other numerical artefacts. Furthermore, numerical results do not 
provide any understanding of the mechanism causing the reduction in dimensionality 
and so we hope the work begun in this paper may be extended to shed some light on this 
question. This is not an unreasonable expectation as the techniques used to compute the 
spectral dimensions of random combs |3j have been extended to allow the computation 
of the spectral dimensions of random trees [18]. Such random trees are closely related 
to two-dimensional CDT via the bijection described in [19] where it was used to bound 
the spectral dimension of the spacetime arising in such models. 

One might ask how closely related a random comb is to an actual model of quantum 
gravity. In models of pure quantum gravity in two dimensions the only observable is the 
spatial volume of the universe. For the case of both CDT and random trees one may 
show that the evolution of the spatial volume with time is generated by a Hamiltonian. A 
crucial difference for the random combs is that the growth of a comb is not a Markovian 
process in the same sense as the growth of a CDT; the vertices at height h from the 
root are not related to those at height h — 1 by a local transfer matrix so a Hamiltonian 
formulation of the evolution of a comb-like universe is not possible. An extension of the 
work described in this paper, using the techniques of [18], to the case of random trees 
would be very interesting as this would constitute an example of dimensional reduction 
for a model which does admit a Hamiltonian formulation. 

We have given in Theorem [10] a reasonably complete classification of what 
behaviours one can find on a continuum comb and indeed it is likely that the cases not 
covered by the theorem do not have a well defined spectral dimension. The behaviour 
for the combs covered by Theorem [10] is fairly rich as there exists a cross-over region 
between the UV and IR behaviours in which a hierarchy of apparent spectral dimensions 
exist. The proof of the bijection between two-dimensional CDT and trees in [19] shows 
that any ensemble of critical Galton- Watson trees is in bijection with a CDT-like theory. 
It is only when the random tree has the uniform measure that we obtain precisely CDT 
on the other side of the bijection. This opens the intriguing possibility of constructing 
a variation of CDT, constructed from an ensemble of random trees with a measure that 
differs from the uniform measure and with a dependence on a length scale. This length 
scale could then be scaled while taking a continuum limit, as we have done in the work 
here. Such a model would likely have intermediate length scales with apparent spectral 
dimensions different from the UV and IR values. We hope to pursue these possibilities 
in future work. 
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Appendix A. Standard generating functions 

We record here a number of standard results for generating functions for random walks 
on combs; the details of their calculation are given in [3]. 
On the empty comb C = oo we have 

Poo(x) = 1 - v^, (A.l) 

p<»(x) -a ^ q+v^-q-^)"- 1 (A2) 
F °° [x} ~ (1 x) (i + V5)» - (i - v^r ' [ } 

G<g(x;n) = (l-x) n ' 2 - (A.3) 

°° V ' ^ V ' (1 + y/x) n ~ (1 - y/x) n K ' 

Note that we can promote n to being a continuous positive semi-definite real variable 
in these expressions; G^(x; n) is then a strictly decreasing function of n and P< n (x) a 
strictly increasing function of n. The finite line segment of length I has 

which it is sometimes convenient to write as 

Pt(x) = v^rtanh (m OC) (x)£) (A. 5) 

where 

WooW = o lo Si /=■ ( A - 6 ) 

2 1 — \/x 



Again, I can be promoted to a continuous positive semi-definite real variable of which 
Pe(x) is a strictly increasing function. The first return probability generating function 
for the comb with teeth of length £ + 1 equally spaced at intervals of n is given by 

P^n{x) = 3 " f {X) - \ [(3 - P,(x) - 2P<"(x)) 2 - 4G£)(x; n) 2 ] 1 .(A.7) 

and is obtained by setting £ = oo in this formula. Pt^ n {x) is a strictly decreasing 

function of t and increasing function of n, viewed as continuous positive semi-definite 
real variables. 

We also need the scaling limits of some of these quantities. They are 

lim a~^G^(a£; a~^v) = <^cosech(z/£ 2") (A. 8) 



and 



lima 2 (i — p i 1 (af)] = --f 2 tanhfpf 2 ) 

4 + 4tanhp^ cothz/^ + tanh 2 p^5 \ (A. 9) 



1 j 
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Appendix B. Bump functions 

A function ip : R — > R is a bump function if ip is smooth and has compact support. We 
will now prove some properties concerning the Mellin transformation of a bump function 
ip, Ai[ip](s), which has support on [a, b] where b > a > 0. 

Lemma 11 The critical strip of the Mellin transform of the nth derivative ofip, ip^, 
is C for all n. 

Proof Recall that the Mellin transform is defined by, A4['0](s) = J °° ip(x)x s ~ 1 dx. Since 
ip has compact support we have, 

M[ip (n) }{s) = / ^ n \x)x s - 1 dx (B.l) 



and since tp is smooth, is bounded on [a, b] by some constant K, so, 

r b 

\M[tp {n) }(s)\ < K x s ~ l dx (B.2) 

J a 

and the RHS is finite for all s G C since b > a > 0. This also shows that Ai[ip( n '](s) is 
holomorphic for all s. 

Lemma 12 Given n e Z+, \M[^]{a + iD)\ < -^.M[|^ (n) |](cx + n) for all s G C. 

Proof Recall from the previous lemma that the critical strip of the Mellin transform of 
"0 and its derivatives coincides with C. We can therefore use the integral representation 
of the Mellin transform to prove statements valid for all sGC. By integration by parts, 

M[ip {n) ){s) = - [ ^ n+l \x)—dx = --M[4> {n+1) ](s + 1) (B.3) 

J a S S 

and therefore, 

MMs) = f^ n \x)x s+n dx = ^i" 1 ^ - M[^)(s + n). (B.4) 

Hence, 

\M M (a + iT) | <^f \^ (n) (*) \*° +n dx = j^M [|^ (n) |] (a + n) . (B.5) 

Given a bump function ~§± e which is always positive, has support [1, 1 ± e] and is scaled 
such that its integral is one, we define the cut-off function to be, 

v± ( x ) = 1 - / m± e (x)dx. (B.6) 

J — oo 

Lemma 13 The critical strip of Ai[i]±](s) is given be Re[s] > 0. The analytic 
continuation of .M-[f]±\( s ) to all s is given by 

M[ri ± ](s) = ~M[*± e ](s + l). (B.7) 

Proof The analytic continuation of Ai [r]±\ (s) may be obtained by applying integration 
by parts to the Mellin transform of n± and recalling by lemma ( ITT1) that A^[ l l/± e ](s) is 
holomorphic everywhere. 
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Starting with 



S±(y) = 5>(*ht(4/) = A*(0) + ^M%±W = MO) + S£\y) (C.l) 
e=o 1=1 



where r)± is the smooth cut-off function introduced in Appendix B and y controls where 



the cut-off occurs we take the Mellin transform of S± (y) with respect to y, 

POO 

M[s£>)(a)= / S<£\y)ir 1 dy 
Jo 



V,(s)M[^ ±e ](s + l)/s (C.2) 



where ^(s) is the Dirichlet series associated to the measure, 



=1 

7«/ 



It is easy to see that the fundamental strip of S± (y) is Re[s] > 0, due to the compact 
support of 7]± on the positive real axis and so the Mellin transform does indeed exist. 
We may now invert the Mellin transform to obtain, 

S±(y) = ,(0) + -L rrtH**^. (C.4) 

This may be computed by rewriting the above as, 

3±{y) = a*(0) + <f y- s V,(s) M ^ ±e]is + 1) ds - R{y) (C.5) 

where the contour C is the rectangle composed of the points {c — zoo, c + zoo, —N + 
zoo, —N — zoo} with N > 0, c such the contour is the right of all poles and 

m = r A( ,) jMi4 (c.6) 

J— N—ioo s 

If has slow growth in the strip —N < Re[z] < then due to lemma ( Tl2l) . which shows 
that [^±e] decays faster than any polynomial as t goes to infinity, the contributions 
from integrating along the contours c + ioo to —N + ioo and from c — ioo to — iV — ioo 
are zero. Furthermore the remainder term R(y) satisfies, 

N p-N+ioo 

\R{y)\ < V — / \V,(s)\\M[* ±e ](s + l)\ds (C.7) 

iv j_ N _ ioo 

and so will only contribute terms of order y N to S±(y). We therefore have, 

S ± (y) = MO) + ^-J y- s V,( s ) M[ *^ ]{s + 1) ds (C.8) 
^ l Jc s 

- R(y) (C.9) 



res 
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where S is the set of positions of the poles of X> M and we have used the fact that 
^V(O) = 1 — MO). Finally, define x±( u ) — Y^a=o f J, {fyv±{l/ u )- By relating this function 
to S±(y) we may write, 

X±{u) = /x(O) + -L / ^( g ) ^g±flfe±l) d5 (c.10) 

-J2(y). (C.ll) 



1 + res 
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